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SUMMARY 


The tenroerature distributions encountered in thin solid wings 
subjected to aerodynamic heating induce thermal stresses that may 
effectively reduce the stiffness of the wing. The effects of this 
reduction, in stiffness were investigated experimentally by rapidly 
heating the edges of a cantilever plate. TOie midplane thermal 
stresses imposed by the nonuniform temperature distribution caused 
the plate to buckle torsionally, increased the deformations of the 
plate under a constant applied torque, and reduced the frequency of 
the first two natural modes of vibration. By using small-deflection 
theory and employing energy methods, the effect of nonuniform heating 
on the plate stiffness was calculated. The theory predicts the gen- 
eral effects of the thermal stresses, but becomes Inadequate as the 
temperature difference Increases and plate deflections become large. 


INTRODUCTION 


One of the structural problems produced by nonuniform heating is 
a change in the effective stiffness of the structure caused by thermal 
stresses. This change in stiffness is not associated with a change in 
the material properties, but depends on the state of stress and may 
occur at stress levels well below those necessary to produce buckling. 
Laboratory demonstrations have shown that the natural frequency of a 
simplified wing structure can be effectively reduced by the nonunifoim 
temperature distribution associated with rapid heating. 

The type of temperature distribution produced by the aerodynamic 
heating of a thin missile wing is shown in figure 1. This figure 
(given in ref. l) shows the variation in temperature across the chord 
of a solid double-wedge airfoil immediately following a Ig acceleration 
to Mach niimber h at 50^000 feet. The temperatures on the surface and 
at the midplsne are shown. Such a temperature distribution produces 
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con^jressive stresses near the leading and trailing edges and tension in 
, the central part of the wing. The effect of these stresses is a reduc- 
tion in the stiffness of the structure that jnay be observed as an 
increased deformation under load or as a reduction in natural frequency. 
No quantitative data exist, however, which relate changes in frequency, 
and hence effective stiffness, to temperature distribution. The project 
to be discussed in this paper was devised to satisfy part of this need. 


TEST PROCEDURE AND RESULTS 


An arrangement was selected that required a sinrole test setup and 
that could be analyzed theoretically with relative ease. The test 
arrangement used- is shown in figure 2. A nonuniform temperature dis- 
tribution resembling that encountered by aerodynamic ally heated wings 
was produced by rapidly heating (with carbon- rod radiators) the longi- 
tudinal edges of a cantilever aluminum- alloy plate l/li- inch thick and 
20 inches square. The temperature distribution along the span of the 
plate was constant except for a slight decrease toward the tip. Fig- 
ure 3 shows how the ter^erature on the edge and at the midchord varied 
during a test. 

In figure 5 the temperature in degrees Fahrenheit of an edge and 
the midchord line is plotted against time in seconds. Heat is applied 
to the plate for about l6 seconds. During this time, the edge tender a- 
ture rises almost linearly to 315° ^ at the peak of the heating cycle. 
When the heating stops, the tenperature drops quickly and slowly levels 
off as the plate cools. The variation in temperature across the chord 
has been shown in figure 2 for a time in the heating cycle (lO seconds), 
at the time of maximum edge temperature ( 16.5 seconds), and during 
cooling (30 secor^s ) . These distributions show that the tenperature 
remains relatively low over the center half of the plate, but rises 
sharply near the heated edges. 

Two investigations were made to determine the effect of the non- 
unlfom temperature distribution on the stiffness. First, the deforma- 
tions of the plate due to thermal stresses were obtained for various 
load conditions . Second, the changes in natural frequency of vibration 
during heating were investigated. 


Deformations 

Figure shows the tip- rotation histories. The angle of rotation 
of the tip in degrees is plotted against time in seconds during heating 
and cooling for no external load and for applied torques of 400 inch- 
pounds in each direction. In each case, the plate deformed by rotating 
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torsionall;^ about the mldchord line. As the plate cooled^ the deforma^ 
tions decreased and the plate returned to its original position. Note 
that the plate underwent an appreciable deformation vithout the appli- 
cation of an external load. This thermal buckling of the plate may be 
a significant factor in the control of missiles having solid fins. She 
plate rotated in the same direction regardless of whether heat was 
applied S3nnmetrlcally or on either edge because of an initial twist. 

The increased deformations, which occurred when a constant torque was 
applied in the seme direction as the initial twist, indicate an approx- 
imate superposition of the deformation induced by the thermal stresses 
on the initial static deflection. When a torque was applied in the 
opposite direction, the plate twisted in the direction of the applied 
torque, but the maximum twist was slightly less. 


Natural Frequencies 

It is difficult to detect the frequency changes of a plate under 
transient heating conditions because of the time required to establish 
resonance. For this reason the present investigation has so far been 
limited to the first bending and first torsion nK3des, inasisich as a 
frequency history could he obtained for these two n»des hy striking 
the plate at regular intervals during the test. 

Figure 5 shows how the frequencies of the first two modes varied 
during the test. Here, frequency of vibration in cycles per second 
has been plotted against time in seconds for the first bending and 
first torsion modes. The first bending frequency decreased from 
19 cycles per second to I 5 cycles per second at the point of maximum 
temiperature gradient. This is a 21-percent reduction in frequency. 

The first torsion frequency begins at 48 cycles per second and drops 
to a minimum of 30 cycles per second - a reduction in natural frequency 
of about 35 percent. As the plate cools, both frequencies return to 
their original values. The small irregularity which occurs at the 
peak of the heating cycle has been observed in all the first- torsion- 
mode tests, but as yet its cause has not been determined. 


THEORETTTAL ANALYSIS AND COMPARISON WITH EKPlRIMliaTAL RESTJIHS 


As a first approach to predicting the changes measured in these 
tests, small-deflection theory has been used. The analytical approach, 
which is a comiblnatlon of the energy methods used to solve buckling and 
vibrational problems, is outlined in the appendix. 
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Deformations 

A comparison of the predicted and actual effect of the temperature 
gradient on the torsional stiffness for the case of no vibration is 
shown in figure 6. Here the tip rotation in degrees has been plotted 
against the edge- to- center temperature difference AT in degrees 
Fahrenheit for an applied torque of 400 inch-pounds. Curves are shown 
for the period of heating (or increasing AT) and cooling (or decreasing 
at). The small-deflection theory gives a reasonable approximation of 
the reduction in stiffness about halfway to the buckling tenperature. 
Above this pointy however, the small-deflection theory indicates that 
the plate deflections would increase more rapidly than actually occurs. 


Natural Frequency 

Figure 7 coapares the results of the calculations with test results 
when there are no external loads and the plate is vibrating in the first 
torsion mode. In this figure the frequency in cycles per second has 
been plotted against temperatiire difference between the edge and the 
midchord line for the periods of increasing and decreasing AT. The 
small-deflection theory again approximates the frequency change about 
halfway to the critical buckling tenperature . Above this value, the 
theory again overestimates the change. This disagreement occurs because 
the distortions have become large and the small-deflection theory is no 
longer valid. When the analysis is extended to include the effects of 
large deflections, more satlsfactoiy agreement is expected at the 
critical AT. 

Figure 8 shows a similar comparison of measured and calculated 
frequency as a function of time. The ratio of actual frequency to the 
initial uni form- temperature frequency is plotted against time in sec- 
onds. The calculated curve indicates that the theoretical critical- 
buckling temperature differential is reached in about I 5 seconds . The 
small-deflection theory predicts that the plate would have lost all its 
stiffness at this point, but this is not the actual case. Since stiff- 
ness is proportional to the square of the frequency, the frequency 
decrease obtained indicates that only about half of the stlfi^ess was lost 
as a result of the Induced thermal stresses. As the plate cools and 
the tenperature difference becomes less than half that required for 
buckling, the theory is again in fair agreement with the test results. 


CONCrUDING REMARKS 


Tests of a cantilever plate have shown that the midplane stresses 
inposed by a nonuniform temperature distribution can effectively reduce 
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the stiffness of the plate. This reduction in stiffness is reflected 
in the increased deformation ’under the action of a constant applied 
torque and also in the reduction of the natural frequency of vibration 
of the first two inodes of the plate. By using small-deflection theory 
and by employing energy methods, the effect of nonuniform heating on 
the plate stiffness was calculated. The theory predicts the general 
effects of the thermal stresses, but is inadequate when the deforma- 
tions become large. An extension of the analysis to account properly 
for large deflections is expected to give more satisfactory results 
near the critical teu^jerature . 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., May 5^ 1-955 • 
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APPENDIX 

SU4ALL-DEFLECTI01T ANALYSIS OF CAMIILSVEE PLAOS SUBJECOSD 
030 HOMJNIFOEM CT'CERATUEE DIS^ERIBUTIOK 


By using small-deflection tiieory and "by making several simplifying 
ass’jmptions, an analysis of the effect of rapid heating on the deforma- 
tions, natural frequencies of vibration, and torsional stiffness of a 
cantilever plate was made. The method used is outlined in this appendix. 
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(GJ)e 

K 

P 

t 


Symbols 

plate length, in. 

arbitrary coefficients of deflection function 
plate half -width, in, 

Et^ 

plate flexural stiffness, 

12(1 - 

Young's mod\ilus, psl 
bending stiffness, lb-in. ^ 

p 

effective bending stiffness, lb-in. 
concentrated load, lb 
critical axial load, lb 

gravitational constant, in./sec^ 

2 

torsional stiffness, lb-in. 

2 

effective torsional stiffness, lb-in. 

arbitrary coefficient of stress function 
distributed load, psi 
plate thickness, in. 



MCA RM L55E20C 


7 


ci^Uy 


ten 5 >erature , ^ 
tenroerature of midchord^ 

difference In teiriperature between edge and mldchord, ^ 
critical teniperature difference, 

potential energy of bending, in-lb 
strain energy due to thermal expansion, In-lb 
potential energy due to external load, in-lb 
energy due to inidplane stresses, in-lb 
strain energy of midplane stresses, in -lb 
kinetic energy of Yibratlon, In-lb 
deflection, in . 

longitudinal coordinate measured from root, in. 
transverse coordinate measured from midehord line. In. 
coefficient of thermal expansion, 

ln.-°F 

arbitrary exponent in tenperature distributlcna 

tip rotation, deg 

tip rotation when ^ = 0, deg 

Poisson's ratio 

specific weight, Ib/cu in. 

longitudinal direct stress, psi 

transverse direct stress, psl 


shear stress, psl 
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j5 stress function 

03 freq.uency, cps 

03q frequency when = 0, cps 


Ten^eratiire Distribution 

The tea^serature is assumed to he constant through the thickness 
of the plate and along the spanj the distribution across the chord is 
represented by the sickle power lav 

T = (1) 

which involves the edge -to -midchord ten^jerature difference AD and 
the power ^ required to describe the measured chordwise temperature 
distribution. The power which varies during a test, determines 
how sharply the temperature rises near the edges of the plate. In 
order to find how ^ varied, a curve of the form given in equation (l) 
was fitted to the chordwise teniperat\are distributions measured at vari- 
ous times dirring a test. At the start of heating ^ is very large, 
but it decreases to about 4 at the peak; edge temperature, continues to 
drop as the plate cools, and reaches a value of about 2 at the end of 
the test (about 55 seconds after the peak; temperature) . The temperature 
distribution given by eq\oation (l) fits the measured distribution rather 
well during heating but becomes increasingly poor as the plate cools. 


Thermal Stresses 


The plate is assumed to be in a state of plane stress and all 
stresses are assimed to be in the elastic range of the material. The 
assumption also has been made that the material properties of the plate 
do not change with temperature. Thermal stresses are given by the 
relationships 


■'X 




y 


~ dx ^ 


( 2 ) 
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where ^ Is the Airy stress funeticai (see ref. 2) assimed to he given 


^ = KEo_^(x 2 - a2^^^y2 > h2^2 

!lMs fmction satisfies the houndary condition that the stresses 
(eqs. 2 ) all vanish on the free houndaries (x = a and y = +h) . 


In order to evaluate the arbitrary constant in the stress function 
(e^- 3 ), the method of mlninium ccaipleinentary energy is used. The expres- 
sion can be written slnply as 



+ 



= 0 


(l|-a) 






+ 2(1 + tj.) 



dx dy 


Oy^dx dy 


m 


The first integral is the usual strain-energy expression and the second 
represents the strain energy due to thermal expanslcaa of the plate, 

(see ref. 3 .) By using eq.mtion (l) for the tenroerature T and equa- 
tions ( 2 ) and ( 3 ) for the stresses, the constant K may he evaluated 

from equation (ha) . When ^ = 2 the result is 


K _ .z2..T22 I 

a 6 (£ + 3)(S + l) 


(5) 


Effect of Tenroerature Gradient on Deflections 
and Natural Frequencies 

To find the effect of teciperature gradient on the deflections and 
natural frequencies, another mlnimimi-energy method is used. Deflections 
are represented by the power series 


00 00 

= 21 21 

m=l n=l 
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( 6 ) 
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vhich satisfies tbe 'boundary conditions (the slope and deflection are 
zero at the root, x = O) . The undetermined coefficients of the deflec- 
tion function may 'be evaluated by minimizing the' change in energy during 
deflection with respect to each of these coefficients. The expression 
can be written as 


3 a 


ij 


- 0 


( 7 a) 


where 


= i/o7.; I&T * {0) * 0 


dx ajrt 


- 1 r f 


a pb 

i 

0 -b 

2,-j. pa pb 




+ 2t- 


dw ^ 

3 QT ^ 


dx dy 


h( 7 b) 


■ ^/o L 






a pb 

I pw dx dy 

0 ^-b 


The first integral, is the energy due to bending, the second is the energy 
imposed by midplane thermal stresses (obtained from eq.s, 2 and 3)/ 'fclie 
third is the kinetic energy of vibration, and the foxurth is the potential 
energy res-ultlng from any external loads applied perpendicular to the 
plane of the plate. (See refs. 4 and 5.) 


Analysis of torsional deflections .- The thermal buckling which 
occuiTed when the longitudinal edges of the cantilever plate were heated 
restilted in torsional deformations which were similar to the deformations 
that took place when a constant torque was applied to the plate tip. Of 
the first two modes of vibration, both of which experienced a reduction 
in natural frequency, the first 'torsion mode vinderwent the largest change. 
Therefore, the remainder of the analysis has been restricted to torsional 
def cmnatlons . Six terms of the deflection function which are antisym- 
metrieal in y have been used. The deflection is then given as 


w = A-[px% + A2^i^x2y3 + 


Ag2x5y + A24x3y3 + 
AjgX^ty + Aj4X*+y5 


( 8 ) 
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When there is no external load or vibration, the solution of eq,uation 
gives the ten^ierature difference necessary to produce thermal buckling 


For a square plate 



this difference is 


(7a) 


ATcr 


-49.98D 

a^tKBa 


(9) 


When the value of K given by equation (5) is used, equation (9) becomes 

_ _ 18 . 56 D (^ + 5)(C + 1 ) 

^=^"a2tEa I ^ 


Since the flexural stiffness D also contains the modulus E, the 
critical temperature difference depends only on the plate geometry, the 
coefficient of thermal ejcpansion, and Poisson's ratio. 


For the case of no external load or heating, the frequency of the 
first torsion mode of 'a square plate = 2^ is found to be 


° a^jt 



( 11 ) 


The frequency given by this relationship is slightly higher than that 
obtained in the experiments. This is due largely to insufficient 
clamping of the plate root. For the purpose of comparing measured and 
calcTilated results, the theoretical plate was assumed to be slightly 
longer than the actual plate so that the measured and calculated fre- 
quencies would be the same at the start of the test. 

If no heat is applied and the plate is not vibrating, the angle of 
twist of the tip (assuming that the tip remains a straight line) 
resulting from a couple of magnitude 2Fb produced by applying a con- 
centrated load F at the tips (a, "irb) of a square plate = 2^ is 

Go = 32.18 ^ (12) 

Approximate interaction equations .- The deflection modes of the 
plate for first-torsion-mode vibration, thermal buckling, and applied 
torques at the tip are all very similar. If these three modes are 
assumed to be identical, the deflection can then be expressed as a 
function with only one arbitrary coefficient, and simple interaction 
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eq,uations can "be obtained. For heating with no external load, the 
ratio of actual frequency to the initial uniform-temperature frequency 
is found to be 



(13) 


For heating with no vibration, the ratio of the angle of twist of the 
tip to the angle at -uniform plate temperature for a constant torque is 


_ 0 _ ^ 1 

®o 1 _ .M- 
^cr 


(lA) 


The frequency of vibration is proportional to the square root of the 
torsional stiffness GJ, and the angle of twist is inversely proportional 
to the stiffness. Hence, in this simple case it is apparent from equa- 
tions (13) and (l 4 ) that 


(GJ)e ^ ^ 
GJ ' 


(15) 


where GJ is the torsional stiffness of the unstressed plate and (GJ)g 
is the effective stiffness of the plate when subjected to thermal stresses. 


It is of interest to note that a relationship similar to equa- 
tion (15) exists for a simple beam subjected to an axial load. The ratio 
of the effective flexural stiffness of the loaded beam to that of the 
unloaded beam is 


(El)e 

El 



(16) 


where F is the axial load (positive for tension, negative for com- 
pression) and Fgj. is the load necessary to produce column buckling. 

The simple expressions given in equations (15) and (l 4 ) very closely 
approximate the results obtained by solving equation (7a) for related 
values of frequency, temperature, and twist. 
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